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Abstract 

A phase transition signature associated with cumulants of the 
largest fragment size distribution has been identified in statistical mul- 
tifragmentation models and examined in analysis of the ALADIN S254 
data on fragmentation of neutron-poor and neutron-rich projectiles. 
Characteristics of the transition point indicated by this signature are 
weakly dependent on the A/Z ratio of the fragmenting spectator source. 
In particular, chemical freeze-out temperatures are estimated within 
the range 5.9 to 6.5 MeV. The experimental results are well repro- 
duced by the SMM model. 

In nuclear multifragmentation studies a special attention has been paid 
to the largest fragment size which is expected to play the role of the or- 
der parameter, and thus provide valuable insight into the phase behavior 
of investigated systems [H [U O HI [5]. In particular, percolation studies 
have shown that a cumulant analysis focused on the skewness, K^, and the 
kurtosis excess, K4^, of the largest fragment size distribution is a valuable 
method to reveal the presence of a phase transition (critical behavior) in 
finite systems [5]. The percolation transition is indicated by = and a 
K4 minimum of about —1 for all system sizes. Events may be sorted ac- 
cording to various measurable quantities that are correlated with the control 
parameter. 

In this work we show that such cumulant features are not restricted to the 
percolation transition but also are observed at a liquid-gas phase transition 
present in statistical multifragmentation models. With this justification the 
cumulant analysis will be applied to the ALADIN S254 experimental data 
obtained with radioactive beams to investigate the presence and isotopic 
dependence of the transition signals. 

It is instructive to look first at the predictions of the thermodynamic 
model which is known to contain a first-order phase transition. The ther- 
modynamic model is a simplified version of the statistical multifragmenta- 
tion model (SMM), allowing to compute the partition function and thus to 
obtain thermodynamic properties of the system |6j. Our calculations were 
performed for the canonical ensemble of noninteracting one-component frag- 
ments at the constant freeze-out volume three times larger than the normal 
nucleus volume. Figure 1 shows the specific heat and the statistical mea- 
sures of the largest fragment size distribution for three systems with 64, 
216, and 1000 nucleons. The broadening of Cy peaks and reduction in peak 
intensities, as well as the decrease of the transition temperatures indicated 
by the Cy maxima, is observed with decreasing system size due to finite-size 
effects. 
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Figure 1: Canonical thermodynamic model predictions. From top to bottom: the 
specific heat, the mean, skewness and kurtosis of the largest fragment mass distri- 
bution as a function of the freeze-out temperature for different system masses. 



At the transition point, the mean A^ax exhibits the fastest decrease, 
which tends with increasing system size to a step discontinuity expected in 
the thermodynamic limit. It is also clearly seen that the transition is as- 
sociated with K3 = and minimum K4, as in the case of the percolation 
transition. Although these cumulant features are similar for both transi- 
tion types, the shapes of the transitional A^ax distributions are distinctly 
different. In the present model the distribution has a bimodal structure 
as expected for the order parameter at a first-order phase transition in the 
canonical ensemble. In percolation the transition is continuous and the dis- 
tribution is single-peaked. 

In the next step we have examined the cumulants of the largest fragment 
charge {Zmax) distributions within the SMM model f7]. Calculations were 
performed with the standard SMM code which includes secondary decays. 
The freeze-out density of one third of the normal nuclear density was as- 
sumed. For a given system (A,Z), events were uniformly generated over a 



wide range of the excitation energy, and then sorted according to the mea- 
surable quantity Zbound used as a control parameter. The simulation results 
have shown that, similarly as in the thermodynamic model, the cumulant 
signal is well observed at a point where Zmax rapidly decreases. The mean 
excitation energy and the mean microcanonical temperature at this Zbound 
point correspond to the flattest part of the caloric curve. The temperature 
derived in this way, referred to as the SMM breakup temperature, T^, is 
shown in Fig. 2(a) as a function of the system charge Z for the A/Z ra- 
tios of 2.17 and 2.49. These ratios are relevant to the neutron-poor and 
neutron-rich systems investigated in the S254 experiment. As one can see, 
the temperature is around 6 MeV and weakly depends on the system size 
and the A/Z ratio. For a comparison. Fig. 2(b) shows the Hartree-Fock 
limiting temperatures [8]. In this case, regarding their absolute magnitude 
which depends on the nuclear potential form used in the calculations [9] , the 
mass/charge dependences are much stronger. The discrimination between 
the two scenarios was one of the motivations for the S254 experiment. 
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Figure 2: SMM breakup temperatures (a) and limiting temperatures from [S] (b). 



The ALADIN experiment S254 was conducted at GSI Darmstadt. 
Beams of ^^''Sn, ^^'^Sn and ^^^La were used to investigate isotopic effects 
in projectile-spectator fragmentation at 600 AMeV. The secondary beams 
with neutron-poor ^"^Sn and ^^^La projectiles contained also some fraction 
of neighbouring isotopes. The mean compositions of the nominal ^"^^Sn 
(^24La) beams were {Z) = 49.7(56.8) and {A/Z) = 2.16(2.19), respectively 

The most prominent result of the experiment is the observation that 
the isotopic dependence of the projectile fragmentation is weak [TT]. The 
mean IMF multiplicity, the mean largest fragment charge, and the double- 



isotope temperatures, ThcU and TbcLi, determined as a function of Zhound, 
are nearly invariant with respect to the projectile A/Z ratio. These can be 
seen in the top panels of Fig. 3 for the former quantities. When comparing 
the Sn and La results, the resemblance is better if {Zmax) and Zf,ound are 
normalized to the projectile charge. Figure 3 allows to examine the pattern 
of 

Zfnax fluctuations by looking at the variance and the K^^ K4 cumulants 
of the Zmax distribution. Also in this case the results are similar for all the 
projectiles. 



N 



N 



N 



20 



E 50- 



107e 


7 X \ Sn 

" °, , , 1 , , , 1 




- Data 
; — SMM 














'....^—r^ \ 























20 



40 



20 



40 



20 



40 



"^bound 



Figure 3: Experimental data versus SMM predictions. From top to bottom: the 
mean IMF multiplicity, the mean, variance, skewness and kurtosis of the Z^ax 
probability distribution as a function of Zbound- 



The signals characteristic of a phase transition (K3 = and minimum 
K4) are observed at Zbound = 28, 29 and 30 for the ^°^Sn, ^^^gn and ^^^La 
projectiles, respectively. It should be noted that the shapes of the investi- 
gated Zmax distributions are highly irregular since the constraints imposed 
by fixed Zbound values are strong in such small systems. Therefore, the 



presence of some distribution features {e.g. bimodality) cannot be trustwor- 
thy concluded. Despite that, the evolution of the cumulant values is quite 
smooth so that the transition signals are clearly identified. Given the signal 
locations on the Zbound axis, we could estimate parameters of the fragment- 
ing sources at the transition as shown in Table 1. The mean charge of the 
fragmenting system, Zq, was estimated on the basis of percolation and SMM 
simulations as 10-20% larger than the Abound value. The values of the iso- 
tope temperature TueLi were raised by 15% with a possible error of ±5%, 
according to the expected range of corrections for secondary decays feeding 
[12] . The obtained temperatures are in an approximate agreement with the 
SMM predictions given in the last column. The SMM breakup temperatures 
were read out from Fig. 2 with the assumption that the A/Z ratio of the 
fragmenting system remains that of the projectile. 
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32.2±2 


5.9±0.3 


5.87±0.05 




29 


33.5±2 


6.5±0.3 


6.18±0.05 


i24La 


30 


34.5±2 


6.2±0.3 


5.82±0.05 



Table 1: Freeze-out characteristics of the fragnienting systems at the transition 
indicated by Zmax fluctuations. 

For more detailed comparisons with the experimental data we have per- 
formed SMM ensemble calculations following the procedure described in 
|13) . The ensemble specifies the probability distribution of the thermalized 
sources defined by A, Z and E* . Parameters of the ensemble were adopted 
from earlier results for ^^^Au fragmentation [13], and then individually ad- 
justed to improve agreements with the experimental IMF multiplicity distri- 
butions. The quality of the obtained agreements can be seen in the top panel 
of Fig. 3. The next panels show the predicted characteristics of the Z^ax 
distributions. The overall agreement with the data can be concluded as sat- 
isfactory, considering the fact that no optimization of the model parameters 
was performed. 

In summary, simulations with the canonical thermodynamic model and 
the SMM model corroborate the percolation suggestion that cumulants such 
as the skewness (K^) and the kurtosis excess (i^4) of the largest fragment 
size distribution are valuable observables in searching for a phase transi- 
tion in multifragmentation. The location of a transition is well indicated by 
K3 = and K4 minimum, even in very small systems. The measurable quan- 
tity Z}jound may be used to sort events as a control parameter. The cumulant 



analysis applied to the ALADIN S254 experimental data has shown that the 
transition signals are observed at nearly the same Zbound values of 28-30 for 
all the projectiles. The corresponding freeze-out temperatures are in the 
vicinity of 6 MeV and the temperature difference between the neutron-rich 
and neutron-poor sources is within a few percent, in agreement with SMM 
predictions. A satisfactory overall description of the experimental Z^ax dis- 
tributions can be obtained with SMM ensemble calculations. The observed 
agreement with the statistical model indicate that in the multifragmentation 
process the accessible phase-space is of dominant importance. 
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